Revealing mitf functions and visualizing allografted tumor metastasis in colorless and immunodeficient Xenopus tropicalis.

Transparent immunodeficient animal models not only enhance in vivo imaging investigations of visceral organ development but also facilitate in vivo tracking of transplanted tumor cells. However, at present, transparent and immunodeficient animal models are confined to zebrafish, presenting substantial challenges for real-time, in vivo imaging studies addressing specific biological inquiries. Here, we employed a mitf-/-/prkdc-/-/il2rg-/- triple-knockout strategy to establish a colorless and immunodeficient amphibian model of Xenopus tropicalis. By disrupting the mitf gene, we observed the loss of melanophores, xanthophores, and granular glands in Xenopus tropicalis. Through the endogenous mitf promoter to drive BRAFV600E expression, we confirmed mitf expression in melanophores, xanthophores and granular glands. Moreover, the reconstruction of the disrupted site effectively reinstated melanophores, xanthophores, and granular glands, further highlighting the crucial role of mitf as a regulator in their development. By crossing mitf-/- frogs with prkdc-/-/il2rg-/- frogs, we generated a mitf-/-/prkdc-/-/il2rg-/- Xenopus tropicalis line, providing a colorless and immunodeficient amphibian model. Utilizing this model, we successfully observed intravital metastases of allotransplanted xanthophoromas and migrations of allotransplanted melanomas. Overall, colorless and immunodeficient Xenopus tropicalis holds great promise as a valuable platform for tumorous and developmental biology research.

Intestinal remodeling during Xenopus metamorphosis as a model for studying thyroid hormone signaling and adult organogenesis.

Intestinal development takes places in two phases, the initial formation of neonatal (mammals)/larval (anurans) intestine and its subsequent maturation into the adult form. This maturation occurs during postembryonic development when plasma thyroid hormone (T3) level peaks. In anurans such as the highly related Xenopus laevis and Xenopus tropicalis, the larval/tadpole intestine is drastically remodeled from a simple tubular structure to a complex, multi-folded adult organ during T3-dependent metamorphosis. This involved complete degeneration of larval epithelium via programmed cell death and de novo formation of adult epithelium, with concurrent maturation of the muscles and connective tissue. Here, we will summarize our current understanding of the underlying molecular mechanisms, with a focus on more recent genetic and genome-wide studies.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs.

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

DNA methylation clocks for clawed frogs reveal evolutionary conservation of epigenetic aging.

To address how conserved DNA methylation-based epigenetic aging is in diverse branches of the tree of life, we generated DNA methylation data from African clawed frogs (Xenopus laevis) and Western clawed frogs (Xenopus tropicalis) and built multiple epigenetic clocks. Dual species clocks were developed that apply to both humans and frogs (human-clawed frog clocks), supporting that epigenetic aging processes are evolutionary conserved outside mammals. Highly conserved positively age-related CpGs are located in neural-developmental genes such as uncx, tfap2d as well as nr4a2 implicated in age-associated disease. We conclude that signatures of epigenetic aging are evolutionary conserved between frogs and mammals and that the associated genes relate to neural processes, altogether opening opportunities to employ Xenopus as a model organism to study aging.

Mutations in the histone methyltransferase Ezh2 drive context-dependent leukemia in Xenopus tropicalis.

CRISPR-mediated simultaneous targeting of candidate tumor suppressor genes in Xenopus tropicalis allows fast functional assessment of co-driver genes for various solid tumors. Genotyping of tumors that emerge in the mosaic mutant animals rapidly exposes the gene mutations under positive selection for tumor establishment. However, applying this simple approach to the blood lineage has not been attempted. Multiple hematologic malignancies have mutations in EZH2, encoding the catalytic subunit of the Polycomb Repressive Complex 2. Interestingly, EZH2 can act as an oncogene or a tumor suppressor, depending on cellular context and disease stage. We show here that mosaic CRISPR/Cas9 mediated ezh2 disruption in the blood lineage resulted in early and penetrant acute myeloid leukemia (AML) induction. While animals were co-targeted with an sgRNA that induces notch1 gain-of-function mutations, sequencing of leukemias revealed positive selection towards biallelic ezh2 mutations regardless of notch1 mutational status. Co-targeting dnm2, recurrently mutated in T/ETP-ALL, induced a switch from myeloid towards acute T-cell leukemia. Both myeloid and T-cell leukemias engrafted in immunocompromised hosts. These data underline the potential of Xenopus tropicalis for modeling human leukemia, where mosaic gene disruption, combined with deep amplicon sequencing of the targeted genomic regions, can rapidly and efficiently expose co-operating driver gene mutations.

Generation of translucent Xenopus tropicalis through triple knockout of pigmentation genes.

Amphibians generally have three types of pigment cells, namely, melanophores (black and brown), xanthophores (yellow and red), and iridophores (iridescent). Single knockout of the tyr, slc2a7, and hps6 genes in Xenopus tropicalis results in the absence of melanophores, xanthophores, and iridophores, respectively. The generation of triple- knockout (3KO) X. tropicalis for these three genes could allow for observation of internal organs without sacrificing the animals, which would be transparent due to the absence of pigments. In this study, we generated 3KO X. tropicalis, which is one of the most widely used model amphibians, through crossing of a slc2a7 single-knockout frog with a tyr and hps6 double-knockout frog, followed by intercrossing of their offspring. The 3KO tadpoles had transparent bodies like the nop mutant and the frogs had translucent bodies. This translucency allowed us to observe the heart, lungs, stomach, liver, and digestive tract through the ventral body skin without surgery. After intravital staining, 3KO X. tropicalis showed much clearer fluorescent signals of mineralized tissues compared with the wild type. These 3KO X. tropicalis provide a useful mutant line for continuous observation of internal organs and fluorescent signals in the body. In particular, such 3KO frogs would revolutionize fluorescence monitoring in transgenic tadpoles and frogs expressing fluorescent proteins.

Identification of tumor-related genes via RNA sequencing of tumor tissues in Xenopus tropicalis.

Cancer treatment is still challenging because the disease is often caused by multiple mutations. Although genomic studies have identified many oncogenes and tumor suppressor genes, gene sets involved in tumorigenesis remain poorly understood. Xenopus, a genus of aquatic frogs, is a useful model to identify gene sets because it can be genetically and experimentally analyzed. Here, we analyzed gene expression in tumor tissues of three individuals in Xenopus tropicalis and identified 55 differentially expressed genes (DEGs). Gene ontology (GO) analysis showed that the upregulated genes in the tumor tissues were enriched in GO terms related to the extracellular matrix and collagen fibril organization. Hierarchical clustering showed that the gene expression patterns of tumor tissues in X. tropicalis were comparable to those of human connective, soft, and subcutaneous tissue-derived cancers. Additionally, pathway analysis revealed that these DEGs were associated with multiple pathways, including the extracellular matrix, collagen fibril organization, MET signaling, and keratan sulfate. We also found that the expression tendency of some DEGs that have not been well analyzed in the cancer field clearly determines the prognosis of human cancer patients. This study provides a remarkable reference for future experimental work on X. tropicalis to identify gene sets involved in human cancer.

Phenotype-genotype relationships in Xenopus sox9 crispants provide insights into campomelic dysplasia and vertebrate jaw evolution.

Since CRISPR-based genome editing technology works effectively in the diploid frog Xenopus tropicalis, a growing number of studies have successfully modeled human genetic diseases in this species. However, most of their targets were limited to non-syndromic diseases that exhibit abnormalities in a small fraction of tissues or organs in the body. This is likely because of the complexity of interpreting the phenotypic variations resulting from somatic mosaic mutations generated in the founder animals (crispants). In this study, we attempted to model the syndromic disease campomelic dysplasia (CD) by generating sox9 crispants in X. tropicalis. The resulting crispants failed to form neural crest cells at neurula stages and exhibited various combinations of jaw, gill, ear, heart, and gut defects at tadpole stages, recapitulating part of the syndromic phenotype of CD patients. Genotyping of the crispants with a variety of allelic series of mutations suggested that the heart and gut defects depend primarily on frame-shift mutations expected to be null, whereas the jaw, gill, and ear defects could be induced not only by such mutations but also by in-frame deletion mutations expected to delete part of the jawed vertebrate-specific domain from the encoded Sox9 protein. These results demonstrate that Xenopus crispants are useful for investigating the phenotype-genotype relationships behind syndromic diseases and examining the tissue-specific role of each functional domain within a single protein, providing novel insights into vertebrate jaw evolution.

Expanding EMC foldopathies: Topogenesis deficits alter the neural crest.

The endoplasmic reticulum (ER) membrane protein complex (EMC) is essential for the insertion of a wide variety of transmembrane proteins into the plasma membrane across cell types. Each EMC is composed of Emc1-7, Emc10, and either Emc8 or Emc9. Recent human genetics studies have implicated variants in EMC genes as the basis for a group of human congenital diseases. The patient phenotypes are varied but appear to affect a subset of tissues more prominently than others. Namely, craniofacial development seems to be commonly affected. We previously developed an array of assays in Xenopus tropicalis to assess the effects of emc1 depletion on the neural crest, craniofacial cartilage, and neuromuscular function. We sought to extend this approach to additional EMC components identified in patients with congenital malformations. Through this approach, we determine that EMC9 and EMC10 are important for neural crest development and the development of craniofacial structures. The phenotypes observed in patients and our Xenopus model phenotypes similar to EMC1 loss of function likely due to a similar mechanism of dysfunction in transmembrane protein topogenesis.

Molecular functions of the double-sided and inverted ubiquitin-interacting motif found in Xenopus tropicalis cryptochrome 6.

Cryptochromes (CRYs) are multifunctional molecules that act as a circadian clock oscillating factor, a blue-light sensor, and a light-driven magnetoreceptor. Cry genes are classified into several groups based on the evolutionary relationships. Cryptochrome 6 gene (Cry6) is present in invertebrates and lower vertebrates such as amphibians and fishes. Here we identified a Cry6 ortholog in Xenopus tropicalis (XtCry6). XtCRY6 retains a conserved long N-terminal extension (termed CRY N-terminal extension; CNE) that is not found in any CRY in the other groups. A structural prediction suggested that CNE contained unique structures; a tetrahelical fold structure topologically related to KaiA/RbsU domain, overlapping nuclear- and nucleolar-localizing signals (NLS/NoLS), and a novel motif (termed DI-UIM) overlapping a double-sided ubiquitin-interacting motif (DUIM) and an inverted ubiquitin-interacting motif (IUIM). Potential activities of the NLS/NoLS and DI-UIM were examined to infer the molecular function of XtCRY6. GFP-NLS/NoLS fusion protein exogenously expressed in HEK293 cells was mostly observed in the nucleolus, while GFP-XtCRY6 was observed in the cytoplasm. A glutathione S-transferase (GST) pull-down assay suggested that the DI-UIM physically interacts with polyubiquitin. Consistently, protein docking simulations implied that XtCRY6 DI-UIM binds two ubiquitin molecules in a relationship of a twofold rotational symmetry with the symmetry axis parallel or perpendicular to the DI-UIM helix. These results strongly suggested that XtCRY6 does not function as a circadian transcriptional repressor and that it might have another function such as photoreceptive molecule regulating light-dependent protein degradation or gene expression through a CNE-mediated interaction with ubiquitinated proteins in the cytoplasm and/or nucleolus.

Genome-wide identification and spatiotemporal expression profiling of zinc finger SWIM domain-containing protein family genes.

The biological function of the novel zinc-finger SWIM domain-containing protein family (ZSWIM) during embryonic development remains elusive. Here, we conducted a genome-wide analysis to explore the evolutionary processes of the ZSWIM gene family members in mice, Xenopus tropicalis, zebrafish, and humans. We identified nine putative ZSWIM genes in the human and mouse genome, eight in the Xenopus genome, and five in the zebrafish genome. Based on multiple sequence alignment, three members, ZSWIM5, ZSWIM6, and ZSWIM8, demonstrated the highest homology across all four species. Using available RNA sequencing (RNA-seq) data, ZSWIM genes were found to be widely expressed across different tissues, with distinct tissue-specific properties. To identify the functions of the ZSWIM protein family during embryogenesis, we examined temporal and spatial expression patterns of zswim family genes in Xenopus embryos. Quantitative real-time polymerase chain reaction (qRT-PCR) revealed that each member had a distinct expression profile. Whole-mount in situ hybridization showed that both zswim1 and zswim3 were maternally expressed genes; zswim5 and zswim6 were expressed throughout embryogenesis and displayed dynamic expression in the brain, eyes, somite, and bronchial arch at the late tailbud stages; zswim7 was detected in the eye area; zswim8 showed a dynamic expression pattern during the tailbud stages, with expression detected in the brain, eyes, and somite; zswim9 was faintly expressed throughout embryonic development. This study provides a foundation for future research to delineate the functions of ZSWIM gene members.

Modification of maternally defined H3K4me3 regulates the inviability of interspecific Xenopus hybrids.

Increasing evidence suggests that interspecific hybridization is crucial to speciation. However, chromatin incompatibility during interspecific hybridization often renders this process. Genomic imbalances such as chromosomal DNA loss and rearrangements leading to infertility have been commonly noted in hybrids. The mechanism underlying reproductive isolation of interspecific hybridization remains elusive. Here, we identified that modification of maternally defined H3K4me3 in Xenopus laevis and Xenopus tropicalis hybrids determines the different fates of the two types of hybrids as te×ls with developmental arrest and viable le×ts. Transcriptomics highlighted that the P53 pathway was overactivated, and the Wnt signaling pathway was suppressed in te×ls hybrids. Moreover, the lack of maternal H3K4me3 in te×ls disturbed the balance of gene expression between the L and S subgenomes in this hybrid. Attenuation of p53 can postpone the arrested development of te×ls. Our study suggests an additional model of reproductive isolation based on modifications of maternally defined H3K4me3.

Age-associated DNA methylation changes in Xenopus frogs.

Age-associated changes in DNA methylation have been characterized across various animals, but not yet in amphibians, which are of particular interest because they include widely studied model organisms. In this study, we present clear evidence that the aquatic vertebrate species Xenopus tropicalis displays patterns of age-associated changes in DNA methylation. We have generated whole-genome bisulfite sequencing (WGBS) profiles from skin samples of nine frogs representing young, mature, and old adults and characterized the gene- and chromosome-scale DNA methylation changes with age. Many of the methylation features and changes we observe are consistent with what is known in mammalian species, suggesting that the mechanism of age-related changes is conserved. Moreover, we selected a few thousand age-associated CpG sites to build an assay based on targeted DNA methylation analysis (TBSseq) to expand our findings in future studies involving larger cohorts of individuals. Preliminary results of a pilot TBSeq experiment recapitulate the findings obtained with WGBS setting the basis for the development of an epigenetic clock assay. The results of this study will allow us to leverage the unique resources available for Xenopus to study how DNA methylation relates to other hallmarks of ageing.

Evolutionary and functional conservation of myeloid differentiation factor 88 (MyD88) in amphibian Xenopus tropicalis.

As a universal adaptor used by most TLR members, the myeloid differentiation factor 88 (MyD88) plays essential roles in TLR-mediated inflammatory response of invertebrate and vertebrate animals, and functional features of MyD88 remain largely unknown in amphibians. In this study, a MyD88 gene named Xt-MyD88 was characterized in the Western clawed frog (Xenopus tropicalis). Xt-MyD88 and MyD88 in other species of vertebrates share similar structural characteristics, genomic structures, and flanking genes, suggesting that MyD88 is structurally conserved in different phyla of vertebrates ranging from fish to mammals. Moreover, Xt-MyD88 was widely expressed in different organs/tissues, and was induced by poly(I:C) in spleen, kidney, and liver. Importantly, overexpression of Xt-MyD88 triggered a marked activation of both NF-κB promoter and interferon-stimulated response elements (ISREs), implying that it may be play important roles in inflammatory responses of amphibians. The research represents the first characterization on the immune functions of amphibian MyD88, and reveals considerable functional conservation of MyD88 in early tetrapods.

Characterization of a novel corticosterone response gene in Xenopus tropicalis tadpole tails.

Corticosteroids are critical for development and for mediating stress responses across diverse vertebrate taxa. Study of frog metamorphosis has made significant breakthroughs in our understanding of corticosteroid signaling during development in non-mammalian vertebrate species. However, lack of adequate corticosterone (CORT) response genes in tadpoles make identification and quantification of CORT responses challenging. Here, we characterized a CORT-response gene frzb (frizzled related protein) previously identified in Xenopus tropicalis tadpole tail skin by an RNA-seq study. We validated the RNA-seq results that CORT and not thyroid hormone induces frzb in the tails using quantitative PCR. Further, maximum frzb expression was achieved by 100-250 nM CORT within 12-24 hours. frzb is not significantly induced in the liver and brain in response to 100 nM CORT. We also found no change in frzb expression across natural metamorphosis when endogenous CORT levels peak. Surprisingly, frzb is only induced by CORT in X. tropicalis tails and not in Xenopus laevis tails. The exact downstream function of increased frzb expression in tails in response to CORT is not known, but the specificity of hormone response and its high mRNA expression levels in the tail render frzb a useful marker of exogenous CORT-response independent of thyroid hormone for exogenous hormone treatments and in-vivo endocrine disruption studies.

Evolutionarily distinct and sperm-specific supersized chromatin loops are marked by Helitron transposons in Xenopus tropicalis.

Transposable elements (TEs) are abundant in metazoan genomes and have multifaceted effects on host fitness. However, the mechanisms underlying the functions of TEs are still not fully understood. Here, we combine Hi-C, ATAC-seq, and ChIP-seq assays to report the existence of multimegabase supersized loop (SSL) clusters in the Xenopus tropicalis sperm. We show that SSL anchors are inaccessible and devoid of the architectural protein CTCF, RNA polymerase II, and modified histones. Nearly all SSL anchors are marked by Helitrons, a class II DNA transposon. Molecular dynamics simulations indicate that SSL clusters are likely formed via a molecular agent-mediated chromatin condensation process. However, only slightly more SSL anchor-associated genes are expressed at late embryo development stages, suggesting that SSL anchors might only function in sperm. Our work shows an evolutionarily distinct and sperm-specific genome structure marked by a subset of Helitrons, whose establishment and function remain to be explored.

Xenbase: key features and resources of the Xenopus model organism knowledgebase.

Xenbase (https://www.xenbase.org/), the Xenopus model organism knowledgebase, is a web-accessible resource that integrates the diverse genomic and biological data from research on the laboratory frogs Xenopus laevis and Xenopus tropicalis. The goal of Xenbase is to accelerate discovery and empower Xenopus research, to enhance the impact of Xenopus research data, and to facilitate the dissemination of these data. Xenbase also enhances the value of Xenopus data through high-quality curation, data integration, providing bioinformatics tools optimized for Xenopus experiments, and linking Xenopus data to human data, and other model organisms. Xenbase also plays an indispensable role in making Xenopus data interoperable and accessible to the broader biomedical community in accordance with FAIR principles. Xenbase provides annotated data updates to organizations such as NCBI, UniProtKB, Ensembl, the Gene Ontology consortium, and most recently, the Alliance of Genomic Resources, a common clearing house for data from humans and model organisms. This article provides a brief overview of key and recently added features of Xenbase. New features include processing of Xenopus high-throughput sequencing data from the NCBI Gene Expression Omnibus; curation of anatomical, physiological, and expression phenotypes with the newly created Xenopus Phenotype Ontology; Xenopus Gene Ontology annotations; new anatomical drawings of the Normal Table of Xenopus development; and integration of the latest Xenopus laevis v10.1 genome annotations. Finally, we highlight areas for future development at Xenbase as we continue to support the Xenopus research community.

CRISPR-Cas9-Based Functional Analysis in Amphibians: Xenopus laevis, Xenopus tropicalis, and Pleurodeles waltl.

Amphibians have made many fundamental contributions to our knowledge, from basic biology to biomedical research on human diseases. Current genome editing tools based on the CRISPR-Cas system enable us to perform gene functional analysis in vivo, even in non-model organisms. We introduce here a highly efficient and easy protocol for gene knockout, which can be used in three different amphibians seamlessly: Xenopus laevis, Xenopus tropicalis, and Pleurodeles waltl. As it utilizes Cas9 ribonucleoprotein complex (RNP) for injection, this cloning-free method enables researchers to obtain founder embryos with a nearly complete knockout phenotype within a week. To evaluate somatic mutation rate and its correlation to the phenotype of a Cas9 RNP-injected embryo (crispant), we also present accurate and cost-effective genotyping methods using pooled amplicon-sequencing and a user-friendly web-based tool.

An Efficient CRISPR-Cas9 Method to Knock Out MiRNA Expression in Xenopus Tropicalis.

In recent years CRISPR-Cas9 knockouts (KO) have become increasingly utilized to study gene function. MicroRNAs (miRNAs) are short noncoding RNAs, 20-25 nucleotides long, which affect gene expression through posttranscriptional repression. As miRNAs are so small and due to the limitations of known PAM sequences, it is difficult to design CRISPR sgRNAs that reproducibly lead to a KO. We have therefore developed a novel approach using two guide RNAs to effectively "drop out" a miRNA. Validation of efficient CRISPR miRNA KO and phenotype analysis included use of q-RT-PCR and Sanger sequencing. To show specificity of the phenotype, we provide a protocol to use miRNA mimics to rescue the KO phenotype.

Divergent subgenome evolution in the allotetraploid frog Xenopus calcaratus.

Allopolyploid genomes are divided into compartments called subgenomes that are derived from lower ploidy ancestors. In African clawed frogs of the subgenus Xenopus (genus Xenopus), allotetraploid species have two subgenomes (L and S) with morphologically distinct homoeologous chromosomes. In allotetraploid species of the sister subgenus Silurana, independently evolved subgenomes also exist, but their cytogenetics has not been investigated in detail. We used a diverse suite of cytogenetic and molecular FISH techniques on an allotetraploid species in Silurana-Xenopus calcaratus-to explore evolutionary dynamics of chromosome morphology and rearrangements. We find that the subgenomes of X. calcaratus have distinctive characteristics, with a more conserved a-subgenome resembling the closely related genome of the diploid species X. tropicalis, and a more rapidly evolving b-subgenome having more pronounced changes in chromosome structure, including diverged heterochromatic blocks, repetitive sequences, and deletion of a nucleolar secondary constriction. Based on these cytogenetic differences, we propose a chromosome nomenclature for X. calcaratus that may apply to other allotetraploids in subgenus Silurana, depending on as yet unresolved details of their evolutionary origins. These findings highlight the potential for large-scale asymmetry in subgenome evolution following allopolyploidization.